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Adaptive Code-Tracking Receiver For Direct-Sequence Code- 
Division Multiple Access (CDMA) Communications Over 
Multipath Fading Channels And Method For Signal Processing 

In A Rake Receiver 
5 Cross-Reference To Related Application 

This application claims priority of European Patent 
Application No. 00300268.0, which was filed on January 14, 
2000. 

Description 
10 Field Of The Invention 

The invention is directed to a method of signal 
processing in a rake receiver and to a direct-sequence code 
division multiple access (DS-CDMA) receiver in general and 
especially to a direct-sequence code division multiple 
15 access (DS-CDMA) receiver for mobile cellular 
communications. 
Prior Art 

Spread-spectrum communications systems have been used 
for many decades, at first because of the excellent anti- 

20 jamming properties, but lately because of several advantages 
in multiple-access scenarios such as cellular 
communications. In a spread-spectrum system, pseudo-noise 
code signals ( PN code) are used to modulate the data signal, 
where the bandwidth of the former is substantially larger 

25 than that of the latter. The PN code can be used to modulate 
the carrier frequency, called frequency-hopped spread 
spectrum (FH) , or it can be multiplied directly with the 
information signal, called direct-sequence spread-spectrum 
(DS) . The resulting signal then has a bandwidth far broader 

30 than that of the original data signal. 

In the following, only DS-spread-spectrum signals and 
systems are addressed in the prior art discussion. One field 
of applications for such systems is e.g. mobile cellular 
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communications. Recent worldwide standardization efforts 
defined DS-spread-spectrum as the transmission method of 
choice for third-generation mobile telecommunication 
systems. Orthogonal PN-codes are assigned to different users 



enabling receivers to detect and decode their users 1 signal 
despite interference from other users. Other impairments in 



much of the transmitted signal energy as possible in order 
to maximize the signal-to-noise ratio. In order to do so, 
the received signal is synchronously correlated with the 
known PN-code to yield the original information signal. In 

is multipath fading channels however, the signal energy is 

dispersed over a certain amount of time due to distinct echo 
paths and scattering. One crucial task of the receiver is 
thus to estimate the channel to improve its performance. In 
realizable digital communication receivers, the principle of 

20 Synchronized detection is employed, see f.i. Heinrich Meyr, 
Marc Moeneclacy and Stefan Fechtel in "Digital Communication 
Receivers: Synchronization, Channel Estimation and Signal 
Processing", John Wiley and Sons, New York, 1998, meaning 
that estimates for channel parameters are computed and used 

25 in the detection process as if they were the true values. 
This generally results in a small loss in performance 
compared to the case where one has perfect knowledge about 
the channel, but that loss becomes negligible if the data 
after channel decoding are compared. If the receiver has 

30 information about the channel profile, one way of gathering 
signal energy is then to assign several correlator branches 
to different echo paths and combine their outputs 
constructively, a structure known as the rake receiver. 



5 



transmitting simultaneously in the same frequency band, 



10 



the mobile environment include the fading multipath channel 
and Doppler frequency shift due to mobile motion. 

Presently employed spread-spectrum receivers gather as 
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The rake receiver has several fingers, one for each 
echo path, and in each finger, the path delay with .respect 
to some reference delay such as a direct or the earliest 
received path, must be estimated and tracked throughout the 
5 transmission. The tracking unit generates estimates for the 
delays of several multipath components by means of an early 
late timing error detector and utilizes these estimates for 
the different delays either to shift the phase of the 
pseudonoise-code (PN-code) , or to compute intermediate 

10 samples of the data signal at the estimated timing instants 
However, tracking of such delays in general suffers 
from interference from the other paths and becomes 
troublesome and difficult, especially when echo paths lie 
close together, meaning that their arrival times at the 

15 receiver differ by values in the order of the chip duration 
T c , being the inverse of the PN-code signal bandwidth. 

Several algorithms have been proposed in the past for 
tracking of the multipath delays. All of these algorithms 
suffer from either high sensitivity to multipath fading and 

20 thus reduced performance in such environments, or 

implementation complexity which is prohibitive for the use 
in mobile terminals where size, production costs and power 
consumption are critical issues. One example of the former 
is the so-called early-late gate timing error detector (EL- 

25 TED) , being a natural solution to a mathematical 

optimization problem for nondispersive, additive white 
Gaussian noise channels (AWGN) (see Heinrich Meyr, Marc 
Moeneclacy and Stefan Fechtel in "Digital Communication 
Receivers: Synchronization, Channel Estimation and Signal 

30 Processing", John Wiley and Sons, New York, 1998, Heinrich 
Meyr and Gerd Ascheid in "Synchronization in Digital 
Communications", vol. 1, John Wiley & Sons, New York, 1990 
and A. J. Viterbi in "CDMA* - Principles of Spread Spectrum 
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Communication", Addison-Wesley Publishing Company, Reading 
(MA) , 1995) . 

A prior art early-late gate timing error detector 
usually is able to track echo path delay as long as the 
5 relative delay between two paths is larger than some 

threshold. However, if the relative delay is smaller than a 
certain threshold value, the early-late gate timing error 
detector is not apt to longer distinguish the two paths. 
Furthermore, tracking of weaker paths (lower power) can be 

10 disturbed by stronger paths -with a relative delay larger 
than the aforementioned threshold. For an example of an 
algorithm with significantly higher implementation 
complexity than the early-late gate timing error detector, 
see R. A. litis, "An EKF-Based Joint Estimator for 

15 Interference, Multipath, and Code Delay in a DS Spread- 
Spectrum Receiver", IEEE Transactions on Communications, 
Vol. 42, No. 2/3/4, February-April 1994. 

Especially in indoor mobile communication scenarios, 
the path delay of the last identifiable echo path is small 

20 compared to outdoor scenarios. All of the signal energy is 
dispersed around the path with the shortest delay and a 
sufficient echo path resolution is critical for adequate 
receiver performance, expressed for instance as the bit- 
error rate (BER) . None of the existing early-late gate 

25 timing error detector algorithms with reasonable complexity 
is able to track adjacent path delays which are separated by 
less than the mentioned threshold. With such a tracking 
system, the required signal-to-noise ratio to achieve a 
certain BER is reduced significantly. 

30 Consequently, it is an object of the invention to 

provide a direct-sequence code division multiple access (DS- 
CDMA) receiver and a method for signal processing in a rake 
receiver having an improved tracking performance, especially 
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in cases where delay differences of arriving signal 
components of multipath signals are in the order of the chip 
duration . 

According to the invention and in order to 
5 significantly reduce the effect of multipath disturbance on 
the timing error detector, several early and late 
correlators are used and their outputs are adaptively 
weighted. 

The timing error detector can be coherent or 
10 noncoherent - in the former case, the adaptive weighting of 
correlators is identical to pre-f iltering the spread- 
spectrum signal with an adaptive finite impulse response 
(FIR) filter, followed by only one correlator. In the 
noncoherent case, two adaptive filters and two correlators 
15 are required. 

The invention is described below in more detail and 
with reference to the appended drawings in view of most 
preferred embodiments . 
Brief Description Of The Drawings 
20 Figure 1 shows the general structure of a digital 

code- 
tracking loop for each finger in a rake receiver 
for DS-spread-spectrum communications; 
Figure 2 shows the conventional coherent early-late 

25 gate 

timing error detector, embedded in a code-tracking 
loop for one finger in a rake receiver; 
Figure 3 shows the conventional noncoherent early-late 

gate 

30 timing error detector, embedded in a code-tracking 

loop for one finger; 
Figure 4 illustrates the concept of adaptively 
weighting 



5 



Baltersee-Fock-S^^Lz-Rittich 1-1-1 

several early and late correlator branches within the 
timing error detector; 

Figure 5 shows a coherent implementation of the new 
adaptive 

5 early-late gate timing error detector, using one 

adaptive FIR filter in the synchronization path 
Figure 6 shows a noncoherent implementation of the new 
adaptive early-late gate timing error detector, 
using two adaptive FIR filters and two 
10 correlators; 

Figure 7 shows the detector characteristic for a 
conventional coherent early-late timing error 
detector; 

Figure 8 . shows one possible detector characteristic 
15 for the 

new modified adaptive early-late timing error 
detector in its coherent version; 
Figure 9 shows a two dimensional plot of the 
characteristics 

20 of outer finite impulse response filter taps; 

Figure 10 shows an error signal S-curve for EL 
extension with 

a zero-crossing at x = ± 2T C ; 
Figure 11 shows an error signal S-curve for EL 
25 extension with 

a- zero-crossing at x = ± T c ; 
Figure 12 shows timing tracking estimates for a 
conventional 

timing error detector based on an indoor scenario 
30 with a spreading factor SF=4, moving speed of 

receiver 10 km/h and a chip duration T c = 
lsec/3, 84*10 6 ;. 
Figure 13 shows timing tracking estimates for an 
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inventive 



simple adaptive 



timing error detector based on an 
with values SF=4 and a moving 
a speed of about 10 km/h and a 



indoor scenario 



recei 



ver having 



chip 



duration T 



lsec/3, 84*10 6 ; 



Figure 14 shows timing tracking estimates for an 
inventive 

MMSE adaptive timing error detector based on an 
indoor scenario with values SF=4 and a moving 
receiver having a speed of about 10 km/h and a 
chip duration T c = lsec/3, 84*10 6 ; 
Figure 15 shows a two dimensional plot of tracking 
performance values for different timing error 
detectors; 

Figure 16 shows the detection path of a rake 

receiver model for a detailed theoretical 
discussion thereof; 

Figure 17 a noncoherent delay locked loop. model for a 
detailed theoretical discussion thereof; 

Figure 18 a coherent delay locked loop model for a 
detailed 



Detailed Description Of The Invention 

The invention is explained in more detail and in view 
of preferred embodiments below where like numerals are 
designating similar elements in different embodiments of the 
invention . 

The invention is part of a tracking device for a 
digital spread-spectrum rake receiver 1 as shown in Figures 
1 and 16. 

The incoming signal 2 is fed through a pulse matched 
filter 3, matched to the transmission pulse. An analog-to- 
digital (A/D) converter 4 generates equidistant samples 5 of 



theoretical discussion thereof. 
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the pulse matched filtered signal and feeds them to a 
digital interpolator/decimator 6, which generates 
intermediate samples 7, 8 of the same signal at the 
estimated timing instants f for each of the N rake fingers 
5 9, 10, 11. 

The interpolated signal is fed to the detection path 7, 
where the correlation with the PN-code 16 and symbol 
detection takes place in each finger 9, 10, 11 as shown in 
more detail Figure 16 where essentially only the detection . 
10 path 7 is shown . 

Furthermore, the signal 8 is fed to several timing 
error detectors 12 which generate error signals 13 that are 
proportional to the remaining path delay error in the signal 
of the respective Rake finger 9, 10, 11. These error signals 
is are lowpass-f iltered in the loop filter to yield estimates 

for the respective path delay, again for each rake finger 9, 
10, 11, which are then fed back to the 

interpolator/decimator 6, closing the tracking loop 15. The 
timing error detector 12 is modified adaptively for each 

20 finger 9, 10, 11, using channel state information such as 

the powers, the signal-to-noise ratios and the delays of all 
paths assigned to the respective rake fingers 9, 10, 11. It 
is noted that instead of interpolating and decimating the 
data signal for each rake finger 9, 10, 11, the PN-code can 

25 just as well be shifted in order to guarantee the correct 
alignment between data samples and PN-code samples. 

Furthermore, due to normally very slowly varying path 
delays, delay tracking of all fingers 9, 10, 11 can be 
implemented in time-sharing using only one dedicated 

30 tracking structure. 

Conventional Timing Error Detectors 

Conventional early-late timing error detectors are 
depicted in Figures 2 and 3 in the coherent and noncoherent 



8 



Baltersee- 



Fock-SI 




z-Rittich 1-1-1 




versions, respectively. The structure downstream of the 
digital interpolator/decimator 6 is replicated for each rake 
finger 9, 10, 11, whereas the interpolator/decimator 6 is 
used jointly. The samples 5 of the incoming signal 2 z(t = 
5 nT s ) with T s being the time interval between two samples are 
interpolated and decimated at the estimated timing delay 
r for each finger 9, 10, 11. Two data streams 7, 8 are 
generated, one is fed to the detection path 7 and the other 
8 is fed to the synchronization path 15. There, the samples 
10 and a shifted version thereof, spaced ±T c /2 with respect to 
the samples in the detection path, are correlated in a 
correlator 17, 18 with the PN-code, yielding early 20 and 
late 19 estimates of the transmitted symbols. 



is late 19 estimates, are subtracted by a subtraction unit 21 
and multiplied by a multiplier 22 with reconstructed 
transmitted symbols 23. 

The reconstructed symbols 23 can be generated by 
complex conjugately multiplying symbol decisions or pilot 

20 symbols with estimates of the channel phase or the channel 
phasor. The real part of the resulting samples, now on 
symbol rate. 1/T, is then the error signal x, 13 and is fed 
through loop filter 14 with lowpass characteristic to yield 
the estimate for the timing delay f . * 

25 In the noncoherent case as shown in Figure 3, the 

square of the absolute value of early and late symbols is 
computed by absolute squaring elements 24, 25 and subtracted 
to yield the error signal x, 13. For references, see also A. 
Wilde in "Extended tracking range delay-locked loop", 

30 Proceedings of the International Conference on 

Communications, Seattle, USA, 1995 and Volker Aue and 
Gerhard P. Fettweis in "A Non-Coherent Tracking Scheme for 
the RAKE Receiver That Can Cope With Unresolvable 



In the coherent case as shown in Figure 2, early 20 and 
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Multipath", Proceedings of the International Conference on 

Communications, Vancouver, Canada, 1999. 

The Inventive Adaptive Timing Error Detector Concept 

The concept of the new adaptive timing error detector 
5 is depicted in Figure 4 for the coherent case. Instead of 
using one early and one late correlator 17, 18, several 
early and late correlators 26, 27, 28 are employed. The 
digital interpolator /decimator 6' again feeds the correct 
samples at rate 1/T C to the detection path 7. The 

10 synchronization path 8 is fed with samples at rate- M/T c . 

The shifted versions thereof in the following branches 
of correlators 26, 27, 28 may be aligned- symmetrically in an 
early-late fashion with respect to the detection path 7, as 
in the conventional case, or alternatively shifted in the 

is causal or noncausal direction. 

In either case the samples in each branch 7, 8 are 
correlated with the PN-code 16 and adaptively weighted such 
as to minimize the effect of other multipaths tracked in the 
rake receiver 1. 

20 For M = 1, Xi = -1, X 2 = 1, this inventive embodiment 

reduces to the conventional timing error detector. The 
signal processing after summation of all branches of 
correlators 26, 27, 28 is identical to the conventional 
timing error detector, both in the coherent and the 

25 noncoherent case. 

A structure with more than two correlators in an early- 
late delay tracking loop was presented in "Extended tracking 
range delay-locked loop", Proceedings of the International 
Conference on Communications, Seattle, USA, 1995. However, 

30 there the weights of the different correlator branches were 
solely chosen so as to maximize the tracking' range of the 
loop. No effort was made to adaptively combat interference 
from adjacent paths and their respective rake fingers. 
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The idea of pre-f iltering the spread-spectrum signal 
was addressed by R. Sampaio-Neto and R. A. Scholtz in 
"Precorrelation Filter Design for Spread-Spectrum Code 
Tracking in Interference", IEEE Journal on Selected Areas in 
5 Communications, Vol. SAC-3, No. 5, September 1985. There, 
however, fixed pre-correlation filters were employed in an 
analog system and no effort was made to specifically and 
adaptively combat multipath interference. 

For the inventive coherent timing error detector, the 
10 ■ adaptive weighting of several correlation branches 

corresponds to the adaptive finite impulse response (FIR) 
filtering model in . Figure 5. 

Due to the linearity of the correlation operation, the 
weighting and summation can be performed prior to 
15 correlating with the PN-code 16, yielding a much simplified 
structure with only one correlator and one adaptive filter 
29, the taps 30 of which being identical to the weighting 
factors 31, 32, 33 in Figure 4. 

For the noncoherent timing error detector 12, an 
20 equivalent model is depicted in Figure 6. The transition 

from a bank of adaptively weighted correlators to a single 
correlator branch cannot be done here due to the nonlinear 
nature of the absolute-square operation in squaring elements 
24, 25. 

25 However, an extension of the conventional non-coherent 

timing error detector is to also include adaptive FIR 
filters 29a, 29b in both early and late correlator 17, 18 
branches and to jointly optimize their weights so as to 
minimize multipath interference. 

3 0 Adaptive FIR Filter Design 

The aim of adaptively filtering the data samples in the 
synchronization path 8 of each finger 9, 10, 11 is to 
minimize effects of other paths of the multipath signal have 
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on the detector characteristic of the path tracked by the 
respective rake finger 9, 10 or 11 of rake receiver 1. 

A typical detector characteristic for the conventional 
coherent timing error detector, also called S-curve, is 
5 depicted in Figure 7. Let's denote the current path and as 
such the current rake finger by its index i, then Figure 7 
shows the expected value of the error signal xi, 13 for an 
AWGN channel, i.e. for a channel including additive white 
Gaussian noise (AWGN) , and a raised-cosine pulse-matched 

10 filter autocorrelation function as *a function of the "x in 
the respective finger 9, 10, 11. 

The timing error detector of rake finger i is 
influenced by adjacent multipaths lying everywhere but in 
the zero-crossings of the S-curve as for a value of zero no 

is additional contribution arises for this path having a delay 
resulting in a zero value. 

Assuming perfect channel knowledge, meaning a perfect 
knowledge of path delays Xi and powers as well as signal-to- 
noise ratios, the detector characteristic for each path of 

20 the received multipath signal 2 can be designed adaptively 

by pre-f iltering the synchronization path samples, either by 
generating nulls at the other path locations, or by jointly 
minimizing an adequate cost function. 

In the first case, an additional path j with a delay 

25 relative to the current path i which coincides with a zero- 
crossing of the S-curve of path i does not contribute to the 
error signal (x, 13) of the actual rake finger with path i 
being associated to the respective actual rake finger 9, 10 
or 11. 

30 For the new adaptive timing error detector, one 

possibility of computing filter coefficients is to perform a 
constrained linear optimization, where the constraints are 
given by the zero-crossing of the S-curve at x = 0, in order 
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to keep the tracker unbiased, and by the normalization of 
the slope of the S-curve in the origin to some constant. 

The task solved by the timing error detector control 34 
shown in Figure 1 is then to find the finite impulse 
5 response ( FIR) of filter 29 X = [X x ... X h ] which jointly 
minimizes the sum 

of the mean-square of the detector characteristic of 
the current path at the locations of all other paths and the 
estimator variance induced by the channel AWGN, being 

10 proportional to the filter energy (in the case of an AWGN 
channel and M = 1). The calculated result depends on the 
length of the FIR filter, L, on the amount and location of 
the other paths and on the Doppler frequency spread. 

Another possibility, with reduced complexity compared 

is to the first inventive embodiment, comprises to extend the 
conventional timing error detector by one weighted 
correlator e.g. as shown in Figure 4 by numerals 28, 33 (or 
FIR filter taps) so as to "shift" the first zero-crossing 
left or right of the origin of the S-curve in Figure 7 to a 

20 desired location where the influence of the contribution of 
the error signal 13 of other paths is reduced. With this 
last mentioned algorithm one adjacent path of the multipath 
signal 2 can be considered in the optimization process, 
resulting in a loss in performance compared to the first 

25 scheme at significantly reduced implementation complexity. 
The magnitude of the additional correlator is a 
function of the desired zero-crossing; it can be stored in 
memory and accessed online, depending on the location of the 
closest (or strongest) adjacent, path. 

30 Between these two inventive embodiments, a tradeoff 

between performance and complexity is possible. 
Description Of Further Preferred Embodiments 

A further preferred inventive embodiment comprises 
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advanced timing tracking algorithms and methods for an UMTS 



algorithms for the channel tap delays in a multipath channel 
environment, in the form of early-late gate delay-locked 
5 loops (EL-DLL) , are error-prone in the case that adjacent 
paths are less than about 1,5. T c apart, T c being the chip 
duration. In that case, the timing estimates of rake fingers 
assigned by a central control unit tend to merge, resulting 
in a loss of lock of both fingers. It can also be shown that 
10 significant performance gains can be expected if the "path 
resolution" can be increased. This is due to the fact that 
the signal energy dispersed on the multipath channel can be 
gathered more efficiently by the rake receiver. 



15 synchronization path of the timing error detector (TED) , the 
expected value of the error signal (S-curve) of the timing 
error detector of each path can be shaped so as to minimize 
interference of adjacent paths, according to some 
optimization criterion. Different optimization methods allow 

20 a tradeoff between performance and implementation 
complexity . 

Detailed Theoretical Approach In View Of Preferred 
Embodiments 



25 must be able to accurately acguire and then track the 

multipath channel parameters for the desired user during the 
entire duration of the transmission. 

The timing tracking algorithms are particularly 
sensitive to a sufficiently large multipath separation, 

30 because if that is not guaranteed, delay estimates in the 
rake receiver belonging to paths lying "close" together in 
the delay domain will interfere with one another, causing 
the estimates to merge. This means a loss-of-lock for both 



rake Receiver. It can be shown that 



conventional tracking 



By adaptively pre-f iltering th 



e signal in . the 



A rake receiver for CDMA communications in the downlink 
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rake fingers, making a reacquisition necessary. Minimum path 
separations were determined to be about 1,5 T c for the 
coherent, data-aided and decision and phasor directed early- 
late timing error detector, and about 1,75 T c for the 
5 noncoherent, non-data-aided early-late timing error 
detector . 

However, these limitations are of a severely reduced 

influence for the inventive systems. 

Especially in indoor scenarios however, showing a low 
io delay-spread, the assignment and subsequent tracking of rake 

fingers lying closer together than the mentioned thresholds 

has shown to yield significant performance gains, i.e. about 

5dB for an indoor scenario with a delay spread of about 2 T c 

and a spreading factor SF of 4. 
is If the acquisition and control unit (path searcher) is 

thus able to assign rake fingers to paths which lie close 

together, an adequate tracking scheme which prevents fast 

loss-of-lock is provided by the invention. 

Transmission Model And Timing Error Characteristic 
20 Signal Model 

The baseband-equivalent transmitted CDMA signal for one 

user can be expressed as 



(1) 



00 



5(0 % Z °r * A mod* ^~ k V 
* = ~°° [ N~ ] 



25 



with 



30 



(2) 



data symbol sent at a time instant n 
spreading factor 

effective spreading sequence length 
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d= [ d 0 .... ^ _j ] . : effective spreading sequence 

g r (t) : root-raised cosine TX pulse 

T c : chip duration, for UTRA, T c =244.14 ns 



5 



* The multipath fading channel is given by 



(3) 



N -I 
P 



Kt;r)= S c At)S(T-r (0) 
1 = 0 



10 



with x being the relative path delay and t the absolute 
time variable, respectively. 

The discrete-time model for the rake receiver is 
depicted in Figures 16. The received signal is given by 



where * denotes convolution and n(t) includes additive white 
20 Gaussian noise (AWGN) and other-user interference. The 

signal is pulse-matched filtered and fed to several rake 
fingers 9, 10, 11, ideally one for each resolvable path of 
the received multipath signal 2. In each finger 9, 10, 11, 
the estimated path delay f 1 is compensated for, followed 
25 by multiplication with the spreading sequence 16 and 

summation over one symbol period. Finally, the decorrelated 
signals from each finger 9, 10, 11 are input to the rake 
combiner or summation unit 37, which maximizes the signal to 
noise ratio (SNR) at the output by multiplying the signals 
30 with the optimal combiner coefficients, being the complex 
conjugate of the estimated instantaneous channel phasors, 
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(4) 



r(t) = s(t) * h{t\ t) + n(t) 
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* 



normalized to the overall channel energy. 
Timing Error Detectors 

The pulse matched-filtered signal z(t) is sampled at 
twice the chiprate and input to a digital interpolator 6, 
5 whose task is to compensate for the estimated timing offset 
r . At its output, the data is demultiplexed: the "even" 
samples are fed to the detection path 7, whereas the "odd" 
samples are input to the tracking path 8 -. we get two data ' 
streams at 1/T C . In the tracking path 8, the data and a 

10 shifted version of thereof are decorrelated with the 

spreading sequence 16. In the noncoherent version, the data 
samples (now on symbol rate 1/T) , are absolute-squared and 
subtracted, yielding the error signal 13. In the coherent 
version, the samples are subtracted prior to correlation, 

is correlated and multiplied with complex conjugate maximum- 
likelihood symbol estimates, created by multiplying symbol 
decisions or pilot symbols with channel phasor estimates. 
The real part of the resulting signal yields the error 
signal 13.' In either case, the error signal is fed through a 

20 loop filter 14 with lowpass characteristic to yield the 
timing estimate . 

Timing Error Detector Error Characteristic 

The expected value of the error signal x, 13 for the 
coherent timing error detector, also called S-Curve, is 

25 depicted in Figure 7 for an AWGN channel and a raised-cosine 
pulse-matched filter autocorrelation function. The S-Curve 
is an open-loop characteristic, meaning that feedback 
effects are not accounted for. The normal point of operation 
is when the estimated delay t lies around the origin. A 

30 small deviation from zero will generate an error signal 

which is used in the feedback loop to drive the delay error 
back to zero. If the channel is nondispersive, meaning that 
there is only one propagation path, then the error signal is 
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a function of that one path only. 

However, if one or several other paths exist in the 
channel model, then they will contribute to the error 
signal, increasing the noise level by a factor proportional 
5 to their tap power and to the square of the value of the S- 
Curve at their location. If the S-Curve is equal to zero at 
those locations, then those additional paths will have no 
influence on the error signal on average. The basic idea of 
minimizing the influence of adjacent paths is thus to 

10 "reshape" the characteristic such that is has zero-crossings 
at the locations of adjacent paths. This, can be achieved by 
pre-f iltering the signal in the synchronization path 8 and 
it will be further outlined in the following discussion. 
Interference Minimization By Adaptive Pre-Fil tearing And 

15 Structure Of The New Timing Error Detector 

The concept of the new adaptive timing error detector 
is depicted in Figure 4 for the coherent case. Instead of 
using one early and one late correlator, several early and 
late correlators are employed. The digital 

20 interpolator/decimator 6 again feeds the correct samples at 
rate 1/T C , to the detection path. The synchronization path 
is fed with samples at rate M/T c . The shifted versions 
thereof in the following correlator branches may be aligned 
symmetrically in an early-late fashion with respect to the 

25 detection path 7, as in the conventional case, or 

alternatively shifted in the causal or noncausal direction. 
In either case the samples in each branch are correlated 
with the PN-code 16 and adaptively weighted such as to 
minimize the effect of other multipaths tracked in the rake 

30 receiver . 

For the coherent timing error detector, the adaptive 
weighting of several correlation branches corresponds to the 
adaptive finite impulse response (FIR) filtering model in 
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Figure 5. Due to the linearity of the correlation operation, 
the weighting and summation can be performed prior to 
correlating with the PN-code 16, yielding a much simplified 
structure with only one correlator 26 and one adaptive 
5 filter 29, the taps of which being identical to the 
weighting factors 31, 32, 33 in Figure 4. For the 
noncoherent timing error detector, an equivalent model is 
depicted in Figure 6. The transition from a bank of 
adaptively weighted correlators to a single FIR/correlator 

10 branch cannot be done here due to the nonlinear nature of 

the absolute-square operation. However, a natural extension 
of the conventional non- 
coherent timing error detector is to also include adaptive 
FIR filters in both early and late correlator branches and 

15 to jointly optimize their weights so as to minimize 
multipath interference . 
Adaptive Filter Design 

The design of the adaptive FIR filters offers a 
tradeoff between performance in terms of lock-in-behavior 

20 and complexity. Under the assumption of perfect channel 

knowledge, meaning path delays and powers as well as signal- 
to-noise ratios, the detector characteristic for each path 
can be designed adaptively by pre-f iltering the 
synchronization path samples, either by generating nulls at 

25 the other path locations, or by jointly minimizing an 
adequate cost function'. 
Constrained Linear Optimization 

For the new adaptive timing error detector, one 
possibility of computing filter coefficients is to perform a 

30 constrained linear optimization, where the constraints are 

given by the zero-crossing of the S-curve at x = 0, in order 
to keep the tracker unbiased, and by the normalization of 
the slope of the S-curve in the origin to some constant. The 
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task is then to find the FIR filter 29 X = [A-i ... Ji L ] 
.which jointly minimizes the sum of the mean-square of the 
detector characteristic of the current path at the locations 
of all other paths and the estimator variance induced by the 
channel AWGN, being proportional to the filter energy (in 
the case of an AWGN channel and M = 1) . The result depends 
on the length of the FIR filter 29, L, on the amount and 
location of the paths and on the Doppler frequency spread. 
For general references on linear optimization, see e.g. D. 
G. Luenberger, "Optimization by vector-space methods", John 
Wiley and Sons, 1969 and R. Fletcher, "Practical methods of 
optimization", John Wiley and Sons, 1987*. Define a cost 
function V(^) ' to be 



15 



(5) 



I^(A) = ||A>1| 2 +|U|| 2 = A T (A T A+ I) 



20 



where the matrix A is (Np x L) and contains the mean-square 
contribution of each path at each FIR filter tap location: 



(6) 



25 



A = 



4.. 

A K 



Ai,j = ci • y • g(Ti-Sj) 



ci is the root-mean square (rms) power of path i, normalized 
to the total power, y is the overall rms signal to noise 
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ratio (SNR) and Sj is the j-th filter tap location, relative 
to the detection path sample. For the example in Figure 10, 
8 = [- 1,5, -0,5, 0,5, 1,5]. x and 8 are normalized to the 
chip duration. 

Due to the fact that the chips are essentially 
uncorrelated after descrambling, ||X|| 2 = S L i=i A*i 2 . I is the ( L 
x L) identity matrix. The aim is now to minimize cost 
function V (X) , subject to the constraints mentioned above: 



io (7) 



min{V(A) : D A = f } = min V(X) 



AczF 



F = D + f + 7V[D] 

is The matrix D is (2 x L) and contains the expression for 

the zero crossing of the timing error detector 
characteristic at x = 0 and its slope at the same location: 




= (8) 

d_ 

dt 



D 2 ,,=!-S(0L (9) 



g(t) is the pulse-matched filter autocorrelation function, 
25 f = [0 c] T describes the two conditions, c being the 

constant slope at the origin. D+ is the pseudoinverse of D 
and N[D] is its nullspace. The idea is now to move from a 
constrained optimization to an unconstrained one. To that 
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end, the first step is to reduce the dimensionality of the 
problem. A singular value decomposition of D yields 



D = PSQ r , 



(10) 



where P is (N x N ) and Q is (L x L) . Sis of the form 



S = 



0 



0 

6 



(11) 



where D is (v x v) , and diagonal (v < min(Np, L) ) . The 
nullspace of D is equal to the range of a matrix H, being 
defined as 



H = [q v +i • ♦ • q n ] 



(12) 



with the q± being the indexed columns of Q. This identity 
yields the new unconstrained optimization problem 



min V(D + f + H 0) 

eeR n " v 



(13) 



Inserting X=D* f+ H 6 into (5) yields 



V(X) = (D + f+ H 0) T {A T A+ lXD + f + H^) 



D + f 



+ 2 f r (D + ) r (A r A+ I) H 6 +0 T H r (A r A+ I) H 6 



D + ff + B r e + - e T ce 

1 2 



(14) 
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m 



This is a quadratic form in 0. Solving for 0 finally yields 

e = - c + b d5) 

5 mit 

B = (2 f T (D + ) r (A r A+ I) H) r 
C = 2 H r (A r A+ I) H 

10 the FIR filter X is then given by X = D + f + H0. An example 
for a resulting S-curve in this case is shown in Figure 8. 
Here, 3 paths have been assumed at relative delays of Ti = 
O, T2 = T c and x 3 = 2T C/ with root mean-square powers of 0 
dB, -10 dB and -20 dB, respectively. 

is It is easily seen that the contributions of the two 

paths adjacent to the one at the origin are nulled out by 
' the zero-crossings of the S-curve at their respective 
locations. It is seen in the magnified region that the 
detector characteristic for path 1 is nearly zero at the 

20 location of the two other paths. It is also noted that the 
lock-in range of the delay tracking loop, being the region 
from t = 0 to the first zero-crossing in either direction, 
is physically constrained by pre-f iltering the 
synchronization path samples. 

2 5 Reduced Complexity Optimization 

A sub-optimal method, however, with significantly 
reduced complexity compared to the first embodiment is to 
extend the conventional timing error detector by one or two 
weighted correlators (or FIR filter taps) so as to "shift" 

30 the first zero-crossing left or right of the origin of the 
S-curve in Figure 7 to a desired location. With this 
algorithm, one adjacent path can be considered in the 
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optimization process, resulting in a loss in performance 
compared to the first scheme at significantly reduced 
implementation complexity. The magnitude of the additional 
correlator is a function of the desired zero-crossing; it 
5 can be stored in memory and accessed online, depending on 
the location of the closest, (or strongest) adjacent path. 

Consider the case where the early-late structure is 
extended by two correlators lying ± 1,5 T c from the 
detection sample, i.e. an extension with M = 1 , . L = 4 
10 according to Figure 4. The "outer" taps of the adaptive FIR 
■ are assumed to be equal in magnitude, but with opposing 
sign. The inner taps are assumed to be ±1, as in the 
conventional timing error detector structure. With this 
structure, a zero-forcing condition for the timing error 

is detector characteristic at some location x can be derived in 
a straightforward manner. The expected value of the error 
signal x is given by 

E[x] = g(T / 2 + r) - g{-T I 2 + r) + X\g\- 3T 12 + r) - g(3T 12 + r 
c c ^ v c 7 c 

20 

(16) 

where g(t) is assumed to be the pulse-matched filter 
autocorrelation function in the AWGN case. Setting to zero 
25 and solving for X yields (for the AWGN case) 



g(-T I 2 + r) - g(T ' I 2 + r) 

Mr) = £ £ (17 

g(-3T I 2 + r) - g(3T I 2 + r) 

c c 



Once the exact structure of the timing error detector 
30 is known (M, L) , this characteristic depends only on the 

desired location of the zero crossing. It is shown in Figure 
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9 for the structure mentioned above. S-curves for optimized 
timing error detectors are shown in Figures 10 and 11, where 
adjacent paths are assumed to be 2 T c and T c away, 
respectively. The taps indicate the placement and size of 
5 the EL-FIR filter taps with respect to the detection path 
sample . 

Tracking Performance 

Simulated path delay estimates for an indoor channel 
model, SF = 4 and a mobile speed of 10 km/h are shown in 

10 Figures 12, 13 and 14 for the conventional timing error 
detector, an adaptive timing error detector with reduced 
complexity early late extension (EL extension) and the 
- above-explained minimum mean squared error (MMSE) adaptive 
timing error detector, designed using the constrained linear 

is optimization and as expressed e.g. in equation 7, 
respectively. 

It is seen that with the conventional timing error 
detector, the first two paths merge fairly quickly and both ' 
rake fingers subsequently track the same resulting path. The 

20 simple early-late extension is already able to prevent a 

loss-of-lock of any finger, only low-rate fading is seen to 
influence the path delay estimates. For the MMSE adaptive 
timing error detector, it is seen that the tracking jitter 
increases for the- first path, a consequence of two facts: on 

25 one hand, the path itself is strong compared to the other 
two, resulting in a large error signal, on the other hand, 
due to the large resulting filter taps and as such an 
increased filter energy, the loop noise is further enhanced 
for this case. 

30 The improved tracking behavior for narrowly spaced rake 

fingers directly translates into a significant performance 
gain, reflected in Figure 15. It is seen that the inventive 
receiver gains at least about 4 dB compared to the 
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conventional timing error detector and that the simple and 
the MMSE optimization perform nearly equally well. An 
exemplary detection loss with respect to the case with 
perfect channel estimation (the leftmost' curve) , resulting 
5 from phasor and timing tracking, amounts to 0,8 dB in this 
particular scenario, equally distributed among the LMS 
Kalman phasor estimation used here and the adaptive timing 
tracking . 

In general, a significant performance improvement 

10 results if the detector characteristics for all rake fingers 
are adaptively modified according to the inventive concepts. ■ 

Using the knowledge of instantaneous multipath delays, 
powers and SNR (or estimates thereof) , filters can be 
designed for each rake finger which minimize multipath 

is interference as described above. On one hand, more rake 

fingers can be assigned within a small delay range, because 
a loss-of-lock will be much less likely due to less or no 
interference from other paths. 

On the other hand, and as a consequence of the first 

20 fact, the effective SNR at the output of the rake receiver 
will be improved significantly, because more signal energy 
can be gathered by the rake receiver. Equivalently , the 
signal to noise ratio (SNR) requirement in order to achieve 
a certain bit error rate (BER) will be reduced, especially 

25 in indoor or pedestrian scenarios. The necessary 

computations can be performed at significantly lower 
processing rates than- the chiprate, due to the fact that 
channel delay dynamics are usually very slow. 

Also, the computation of the adaptive filter 

30 coefficients in the case of constrained linear optimization 
can be divided in one part depending solely on the structure 
of the timing error detector, i.e. the oversampling factor M 
and the filter length L, and one part depending on the 
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actual instantaneous channel. The former can be computed 
offline once the structure is known, whereas the latter 
should be computed in an online fashion. 

A tradeoff between computation complexity and 
adaptivity or adapting ability is possible if e.g. the 
filter length is kept variable. Also, in contrast to 
existing solutions, no estimates of instantaneous channel 
multipath amplitudes and phases are required according to 
the invention . 



